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PROJECT  SUMMARY 

Reducing  the  drag  of  cannon  launched  artillery  shells  is  a 
potential;  method  of  increasing  the"' range,  of  the  shell.  The' base 
drag  of  the  shell  is  directly  related  to  the  base  pressure  that 
exists  on  the  base  region  area.  A  potential  method  to  reduce  the 
base  drag  is  to, "employ  a  solid  propellant  gas  generator  located 
in  the  shell  aft  end  that  injects  gas  into  the  base  region. The 
mass  injected  can  be  distributed  in  a  number  of  ways.  It  can  be 
introduced  through  the  center  of  the  projectile,  near  the  edge  of 
the  projectile  or.  a  combination  of  these  techniques.  .  •- '  ■  -  .  ■ 

The  Phase  I  SBIR  effort  concentrated  on  analyzing  ^-the 
effects  of  central  injection  or  edge  injection  and r its-  effect  on 


the  base  drag 
The  injectan# 
and  hydrogen, 
for  the  air 


,  The  in jectantj wes“  injected  at  subsonic  velocities* 
was  considered  to  be  either  air  or  a  mixture  of  air 
The  effect  of  injectant  temperature  was  considered 
injection  case.  The  air  hydrogen  mixture  was 


considered  to  be  both  non  reacting  and  also  >  reacting. 


f  lame 


sheet  combustion  model  was  -employed  for  the  reacting  flow  case. 
The  free  stream  Mach  numbers  considered  were  1.4,  1.8,  and  2.2. 
The  "'effect  of  spin  on  thebase  drag  was  also-' investigated  by 
including  a  spin  term  in  the  axisymmetric  flow  equations. 

The  results  obtained  indicate  that  subsonic  base  injection 
can  be  beneficial  in  reducing  base  drag.  The  •  use  of  edge 
injection  gives  higher  values  of  base  drag  reduction  (by 
approximately  20%)  than  center  injection.  Use  of  higher 
temperature  injectant  gas  also  gives  larger  values  of  base  drag 
reduction.  It  appears  thatjmall  injection  amounts  of  burning  gas 
in  the  base  region, Is  more  beneficial  than  large  amounts,  which 
might  blow  off  th^  base  bubble. 

Details  of  the  optimum  injection  scheme,  however,  still 
remain  to  be  determined.  In  particular,  effects  of  flow  channels 
located  in  the  base  region  that  direct  the  flow  either  radially 
inward  or  outward  may  increase  or  possibly  decrease  the  base 
pressure  and  therefore  alter  the  range  of  the  artillery  shell. 
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ABSTRACT 

Reducing  the  drag  of  cannon  launched  artillery  shells  is  a 
potential  method  of  increasing  the  range  of  the  shell.  The  base 
drag  of  the  shell  is  directly  related  to  the  base  pressure  that 
exists  on  the  base  region  area.  A  potential  method  to  reduce  the 
base  drag  Is  to  employ  a  solid  propellant  gas  generator  located 
in  the  shell  aft  end  that  injects  gas  into  the  base  region  The 
mass  injected  can  be  distributed  in  a  number  of  ways.  It  can  be 
introduced  through  the  center  of  the  projectile,  near  the  edge  of 
the  projectile  or  a  combination  of  these  techniques. 

The  Phase  I  SBIR  effort  concentrated  on  analyzing  the 
effects  of  central  Injection  or  edge  injection  and  its  effect  on 
the  base  drag.  The  injectant  was  injected  at  subsonic  velocities. 
The  injectant  was  considered  to  be  either  air  or  a  mixture  of  air 
and  hydrogen.  The  effect  of  injectant  temperature  was  considered 
for  the  air  injection  case.  The  air  hydrogen  mixture  was 
considered  to  be  both  non  reacting  and  also  reacting,  A  flame 
sheet  combustion  model  was  employed  for  the  reacting  flow  case. 
The  free  stream  Mach  numbers  considered  were  1.4,  1.8,  and  2.2. 
The  effect  of  spin  on  the  base  drag  was  also  investigated  by 
including  a  spin  term  in  the  axisymmetric  flow  equations. 
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INTRODUCTION 


Reducing  the  drag  of  a  cannon  launched  artillery  shell  is  a 
potential  method  of  increasing  the  range  of  the  shell.  The  base 
drag  of  the  shell  is  directly  related  to  the  base  pressure  that 
exists  on  the  base  region  area.  A  generic  artillery  shell  is 
illustrated  in  Figure  la.  By  contouring  or  boat  tailing  the  back 
end  of  the  shell  as  shown  in  Figure  lb,  drag  is  reduced  since  the 
pressure  that  occurs  on  the  forward  facing  base  region  area  is 
not  as  low  as  the  base  pressure  and  causes  a  positive  axial 
pressure  force.  An  alternative  approach  is  to  employ  a  solid 
propellant  gas  generator  that  injects  gas  into  the  base  region  as 
illustrated  in  Figure  lc.  The  mass  injected  can  be  distributed  in 
a  number  of  ways.  It  can  be  introduced  through  the  center  of  the 
projectile,  near  the  edge  of  the  projectile  or  a  combination  of 
these  techniques. 

The  effort  to  be  reported  on  was  purely  an  analytical  study 
and  concentrated  on  central  and  edge  mass  injection.  The 
equations  employed  in  the  study  were  the  axisymmetric  Navier 
Stokes  equations.  In  addition  to  the  theoretical  studies 
performed  some  code  modifications  were  required  and  some  relevant 
software  for  pre  and  post  processing  was  developed.  Prior  to 
discussing  the  theoretical  results  the  pre  and  post  processor 
software  will  be  discussed  then  the  code  modifications  and 
finally  the  results  obtained. 


GRID  GENERATOR  PROGRAM 

An  IBM  PC/AT  was  used  to  develope  the  pre  and  post  processor 
computer  programs.  The  programs  were  written  in  Fortran  and  used 
the  IBM  professional  Fortran  compiler.  The  IBM  AT  was  also  used 

to  communicate  with  the  CRAY  X/MP  at  the  Ballistic  Research 

Laboratory  using  a  1200  baud  modem  and  a  terminal  emulator 

software  program  call  PC  Plot. 

The  transformation  which  generates  the  physical  grid  was 

modified  to  allow  for  grid  point  clustering  in  the  horizontal  and 
vertical  direction.  Several  grid  clustering  transf ormat ions  were 
incorporated  into  the  program.  They  consisted  of  an  axial  grid 
clustering  at  one  x  location  and  a  choice  of  three  vertical 
clustering  techniques.  These  are 

1.  a  cosine  clustering 

2.  clustering  near  the  centerline 

3.  clustering  at  one  vertical  y  position 

A  pre-processor  graphics  program  for  grid  visualization  was 
constructed.  This  program  was  used  to  determine  if  suffucient 
grid  resolution  is  achieved  for  the  particular  case  of  interest. 
The  original  mode  of  operation  for  this  grid  program  was  to 
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submit  the  NASJIN  code  to  the  main  frame  as  a  batch  job 
for  an  iteration  of  one  time  step.  The  grid  generated  was  then 
downloaded  to  the  IBM-PC  and  plotted  using  the  grid  program.  This 
method  of  operation  became  cumbersome  when  modifications  to  the 
grid  such  as  clustering  parameters  were  used  to  change  the  grid 
resolutuon.  This  mode  of  operation  was  then  replaced  by  using  the 
IBM  AT  to  perform  the  grid  generation.  The  same  grid  generator 
that  was  in  the  NASJIN  code  was  implemented  on  the  IBM  AT.  After 
selecting  the  appropiate  x  and  y  clustering  values,  these  values 
were  then  used  in  the  input  to  the  NASJIN  code.  This  eliminated 
the  submission  of  a  number  of  one  iteration  cases  with  different 
clustering  values  and  the  downloading  of  the  resulting  grids. 
Also  this  approach  minimizes  computer  charges  incurred  on  the 
main  frame. 

Another  reason  for  doing  the  pre  and  post  processing 
graphics  on  the  IBM  PC  is  to  make  the  graphics  software  mainframe 
independent.  Thus  when  the  NASJIN  code  is  ported  to  another 
mainframe  the  graphics  will  not  have  to  be  modified  for  that 
mainframe . 

GRAPHICS  POST  PROCESSOR  PROGRAM 

To  display  resulting  flow  field  quantities  a  plotting 
package  was  created.  This  allows  visual  examination  of  the 
resulting  flow  field  predicted  by  the  NASJIN  code.  The  plots  that 
are  possible  with  this  program  are  contour  plots,  velocity  vector 
plots  and  flow  field  profile  plots  on  constant  x  or  y  grid  lines. 
The  program  also  has  the  ability  to  allow  selecting  windows  to 
view  the  flow  field.  This  option  permits  enlarging  sections  of 
the  overall  flow  field  to  examine  details  of  a  particular 
flowfield  region.  These  changes  have  not  yet  been  checked  out 
with  the  plot  output  from  the  NASJIN  code.  The  use  of  this 
graphics  program  would  allow  determining  the  extent  and  location 
of  flow  characteristics  such  as  separated  flow  and  shock  waves. 
This  graphics  program  uses  a  software  package  that  runs  on  the 
IBM  AT  called  Multi-Halo  and  is  a  commercial  graphics  package. 

The  method  of  operation  for  the  graphics  program  was  to 
submit  a  case  or  a  number  of  cases  to  the  CRAY.  After  these  cases 
ran  the  plot  file  generated  by  thes  cases  would  be  downloaded  to 
the  IBM  AT  and  the  plots  generated  on  the  screen.  Hard  copies  of 
these  plots  would  then  be  obtained  by  using  a  pen  plotter.  In 
attempting  to  do  this  two  problems  areas  surfaced.  The  first  was 
the  size  of  the  plot  file.  Although  the  plot  file  is  not  large, 
to  download  the  file  would  take  on  the  order  of  an  hour  on  a  1200 
baud  line.  The  other  difficulty  encountered  was  the  noise  on  the 
phone  line.  During  attempting  to  download  a  plot  file,  noise 
would  come  across  the  screen.  It  is  not  known  if  this  noise 
interferred  with  the  data  that  was  transmitted  or  was  only 
related  to  the  display  on  the  screen.  Another  possibility  is 
that  the  PC-PLOT  software  was  inadequate  in  filtering  of  the 
noise. 
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ADAPTIVE  GRID 


The  incorporation  of  an  adaptive  grid  capabilty  into  the 
NASJIN  code  (in  fact  any  CFD  code)  would  enhance  the  codes 
ability  to  resolve  regions  of  high  gradients.  Also  it  could 
reduce  computer  cpu  time  by  allowing  use  of  fewer  grid  points, 
since  sufficient  grid  points  would  be  relocated  to  high  gradient 
regions  from  low  gradient  regions.  In  order  to  accomplish  this  a 
copy  of  the  adaptive  grid  program  was  obtained  from  NASA  AMES 
Research  Center  (1). 

The  adaptive  grid  program  obtained  utilizes  an  existing  flow 
field  solution  to  relocate  grid  points  according  to  local  flow 
gradients.  These  flow  gradients  are  utilized  in  a  tens  ion/tors  ion 
spring  analogy  to  redistribute  grid  points.  Spring  tension  and 
torsion  coefficients  are  additional  inputs  required  by  the  code. 
An  iterative  (human  in  the  loop)  process  is  required  to  determine 
the  correct  combination  and  magnitude  of  these  coefficients  to 
achieve  an  acceptable  solution  from  the  adaptive  grid  code.  The 
end  objective  for  the  adaptive  grid  code  was  to  have  the  ability 
to  run  the  adaptive  grid  code  as  a  subroutine  in  the  NASOIN  code. 
The  original  grid  and  metrics  are  provided  to  the  NASJIN  code, 
and  these  metrics  are  recomputed  during  the  solution  process  at 
specified  iteration  values. 

There  are  two  versions  of  the  NASA  Ames  two  dimensional 
adaptive  grid  code.  The  first  and  older  version  did  not  have  the 
ability  to  perform  grid  adaption  unles  it  was  used  in  a  manual 
(person  in  the  loop)  mode.  The  second  version  had  a  limited 
capability  to  perfrom  self  adaption  by  adjusting  the  required 
spring  and  torsion  constants  in  only  one  direction.  Incorporat ing 
either  adaptive  grid  version  into  the  NASOIN  code  required  that 
the  NASOIN  code  be  stopped  during  the  solution  proceedure,  modify 
the  existing  grid  to  redistribute  points  in  regions  where  flow 
gradients  are  greatest  and  then  restarting  the  code  with  this  new 
grid  and  flowfield.  An  interface  routine  was  required  between  the 
adaptive  grid  subroutine  and  the  NASOIN  code.  This  necessitated 
the  development  of  an  additional  subroutine  which  would  evaluate 
the  transformation  metrics  for  the  modified  grid. 

Alteration  of  the  adaptive  grid  routine  was  necessary  to 
make  allowances  for  the  additional  number  of  dependent  variables 
computed  by  the  NASJIN  code  (the  specie  mass  fraction).  Coding 
was  added  so  that  each  of  the  flow  field  variables  are  determined 
by  interpolation  from  the  modified  grid. 

Prior  to  starting  this  effort,  the  inclusion  of  the  adaptive 
grid  capability  into  the  NASOIN  code  was  believed  to  be  a  simple 
matter  of  employing  the  NASA  Ames  adaptive  grid  code  as  an 
additional  subroutine.  The  function  of  this  subroutine  would  be 
to  automatically  adjust  the  required  input  values.  After 
initiating  this  effort,  it  was  realized  that  these  constants  were 
highly  problem  dependent.  If  incorrect  values  were  used  the 
resulting  grid  lines  would  cross.  What  was  believed  to  be  a 


simple  matter  of  addition  of  a  few  adaptive  grid  related 
subroutines  became  in  fact  a  separate  developement  effort.  The 
effort  required  deteriming  the  appropiate  values  for  the  spring 
and  torsion  constants  for  the  problem  of  interest  and  then 
modifying  the  adaptive  grid  code  to  adjust  these  values  within 
specified  ranges  during  the  solution  of  the  flow  field.  It  was 
not  possible  to  accomplish  this  during  the  Phase  I  effort. 

A  partial  success  for  the  adaptive  grid  code  was,  however, 
was  obtained.  The  adaptive  grid  code  was  successfully  applied  to 
the  flow  over  a  bluff  base  having  a  centered  jet.  The  freestream 
conditions  utilized  were  Mach=1.5,  14.7  psi  at  520  degrees 
Rankine.  The  jet  conditions  were  Mach  1.85,  29.4  psi  (exit 
pressure)  at  630  degrees  Rankine,  the  injectant  was  air.  The  flow 
solution  was  obtained  at  500  time  steps  and  a  restart  tape  was 
written.  This  restart  file  along  with  other  required  inputs 
served  as  input  data  to  the  adaptive  grid  routine  and  the  flow 
was  then  solved  for  another  500  time  steps.  A  plot  of  the  grid  at 
time  steps  of  500  and  1000  Iterations  is  presented  in  Figures  2 
and  3 . 


FIGURE  2 


ADAPTED  GRID  AT  500  ITERATIONS 
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FIGURE  3 

ADAPTED  GRID  AT  1000  ITERATIONS 
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CENTERLINE 


SPIN  MODIFICATIONS 

Analizing  the  effects  of  spin  on  the  base  flow  was  achieved 
by  incorporat ing  additional  terms  into  the  two  dimensional 
axisymmetric  flow  equations.  These  terms  had  the  effect  of 
imposing  an  angular  velocity  everywhere  throughout  the  flow 
field.  This  is  an  approximation  and  represented  an  alternative 
approach  to  solving  the  full  set  of  three  dimensional  flow 
equations.  If  the  three  dimensional  flow  equations  were  solved  a 
constant  angular  velocity  boundary  condition  would  have  been 
imposed  only  on  the  body  surface.  Thus  the  flow  field  model 
employed  for  incorporat ing  spin  effects  is  oversimplified  and 
does  not  represent  the  true  effects  of  spin.  It  does,  however, 


represent  a  worst  case  were  the  spin  effects  would  be 
exaggerated.  A  more  realistic  approach  would  have  been  to  include 
a  0  momentum  equation  and  solve  for  the  velocity  component  in 
this  direction.  The  three  dimensional  Navier  Stokes  equations  in 
a  cylindrical  co-ordinate  system  were  simplified  by  neglecting 
any  flow  variation  in  the  c ircumf erent ia 1  direction  and  assuming 
a  constant  spin  velocity  that  is  independent  of  axial  location. 
This  simplification  leads  to  the  following  equations: 


continuity: 


dp  3 pu  dpu  pv  _ 

dl  +  dx  d y  y 


x-momentum: 


3/,  2  /  3  u  .  3u  i  u  \  _  3u\  .  d  (  (du  3iA  ,  pUu  a  fdu  dv\ 

■SV*1  +p  + in  3i+7  +  +3?H-Ka?  +  3?J+  7  y  U  +  W  ~ 


'V 


y-momentun,  g+  A(?„„_(1(|  +  |))  +  A^.  +J>- 1„(  |-|-  *)-2„  |)  + 


9-momentum 


where : 


^(pmv)  +  |;(pw)+  ^pw=R0 


„  3  /  3w\  3  /  w\  ,  2  /  3w  w\ 

Re-a?l"3?)~3?l''~7)  +7^^r-  *■  7J 


The  next  step  is  to  simplify  the  0  direction  equation  and  then 
obtain  a  suitable  functional  form  for  the  angular  velocity. 

Rearranging  the  terms  on  the  0  equation  there  is  obtained: 


»  ^  +  -22)  +  p»  £  +  m  £  +  -  K 

\  bx  3 y  y  /  3y  3*  v  € 


the  first  term  is  zero  since  it  is  the  steady  state  continuity 
equation.  This  leaves  the  equation: 


(3  w  hA  _ 

by  y  )~ 


Assuming  that  the  w  velocity  component  is  of  the  form: 


w  =  Ay" 

and  substituting  this  into  the  0  equation: 


rho  v  (  An  y"-1  +  Ay"-1)  =  R„ 

in  order  for  the  term  on  the  left  side  of  the  equation  to  be 
zero  n  *  -1  and  A  must  be  a  constant.  The  constant  A  can  be 
determined  using  the  velocity  at  the  body  surface  Wt> . 

wt,  =  fl  yt,  =  A/y*, 


m  XVWTm  JTU  ITJ 


where  ft  is  the  angular  spin  rate.  This  gives  for  the  value  of  A: 


A  =  ft  y„a 


The  end  result  for  the  circumf erent ia 1  velocity  equation  is  then: 


ft  yt=3/y 


This  equation  should  also  be  a  solution  of  the  viscous  RHS  of  the 
theta  equation.  Substituting  this  into  the  RHS  does  in  fact 
indicate  that  it  satisfies  the  equation.  On  the  centerline  this 
equation  is  indeterminate,  therefore  when  the  vertical 
co-ordinate  is  less  than  the  body  radius  it  was  assumed  that  the 
fluid  undergoes  solid  body  rotation.  The  form  for  w  used  in  this 
region  is  : 

w  =  0  r : 

Using  these  forms  for  the  w  velocity,  at  the  body  surface  the  w 
velocity  is  continuous.  The  equations  used  in  the  NASJIN  code 
employed  this  w  velocity  equation  to  asses  the  effects  of  spin  on 
base  pressure. 


EXPERIMENTAL  VERIFICATION 

The  experimental  verification  of  results  predicted  by  the 
NASJIN  code  would  serve  to  establish  a  confidence  level  for  the 
code.  To  accomplish  this  an  experimental  program  that 
investigates  the  effects  of  base  injection  location  (i.e  center 
or  edge  injection)  and  the  effects  of  combustion  would  be 
required.  A  survey  was  made  to  determine  the  availability  of  wind 
tunnel  facilities  where  injection  tests  with  combustion  could  be 
performed.  The  objective  of  the  experimental  program  is  to  gather 
detailed  flow  field  data  necessary  to  validate  the  code.  The 
experiments  would  be  designed  to  supplement  deficiencies  in 
available  test  data.  The  primary  goal  is  to  define  the  pertinent 
flow  and  geometry  parameters  which  minimize  the  projectile  base 
drag.  A  literature  survey  was  performed  to  determine  the  extent 
of  existing  data  relevant  to  the  problem  of  interest.  The  survey 
showed  that  there  is  substantial  experimental  scatter  and  lack  of 
agreement  between  prior  theoretical  models.  Also  the  influence  of 
several  important  flow  parameters  such  as  discrete  injection,  and 
model  spin  is  either  totally  lacking  or  does  not  cover  a 
sufficiently  wide  range  of  interest.  Available  test  facilities 
for  full  scale  experiments  were  examined  on  the  basis  of 
suitability,  availability  and  cost.  A  list  summarizing  the 
results  of  the  available  wind  tunnel  facilities  is  given  in  Table 


WWScwP 


appropiate  faciltly  to  perform  the  tests  is  that  at  General 
Applied  Science  Laboratories. 


Table  l 


Test  Facilities  Contacted 


Facility 

/Location 

Test 

Section 

Flow 

Conditions 

Comments 

Cost/ 

Basis  . 

General  Applied 
Science,  N.Y. 

8  “  x  1 0  " 

M  «  2.7 

Combustion 

Testing 

$  9000/wk 

Naval  Surface 
Weapons  Center 
lid. 

16 "xl 8" 

. 3  <M<  5 .0 

Half  Scale 

No  combustion 

-- 

NASA  LRC 

Va. 

4  '  x4  ' 

1 .47<M<2 .86 

No  combustion 

Grumman 

Aircraft,  N.Y. 

15" 

M  *  2.2 

Half  Scale 

Wright  Aero. 

Lab .  ,  Ohio 

8 . 2 " x  8" 

M  »  3.0 

No  combustion 

2000/day 

LTV  Texas 

4  'x4  ' 

No  combustion 

1375/hr 

Rockwel 1 

Internet iona  1 
Ca. 

7  '  x7  ' 

M  =  2.0 

Some 

combustion 

2800/hr 

National 

Research 

Council,  Can. 

5  '  x5  ' 

. 1 <M<  4.25 

Combustion 

12500/day 

NASA  LeRC 

1  'xl  ' 

1 .6<M<5 .0 

— 

— 

NASA  ARC 

6  '  x6  ' 

1  .  5  <li<  2 . 5 

No  combustion 

— 

AEOC 

APTU 

16 'xl6  ' 

0 . 0  <M<  4 . 5 

Combustion 

— 

New  York  Univ. 

8  "  x  1 0  " 

M  «  2 .7 

No  combustion 

2000/day 

S! 


PROBLEM  GEOMETRY 


The  geometry  employed  for  the  center  injection  cases  is 
illustrated  in  Figure  4a  and  the  geometry  employed  for  the  edge 
injection  cases  Is  illustrated  in  Figure  4b.  The  number  of  grid 
points  used  for  each  of  these  configurations  is  listed  in  the 
figure.  Typically  100  points  were  used  in  the  axial  direction  and 
50  grid  points  in  the  vertical  direction.  A  cartesian  grid  was 
specified  with  suitable  clustering  transformat  ions  employed  to 
resolve  flow  details  in  the  vicinity  of  the  corner  and  injector 
areas . 


55  Grid  Points 


UPPER  BOUNDARY 


CENTER  INJECTION  GEOMETRY 


INFLOW  PLANE 


OUTFLOW  BOUNDARY 


100  Grid  Points 


39  Grid  Points 


21  Grid  Points 


INJECTANT  NOZZLE, 


x  -  inches 


5 1  52 


CENTER  INJECTION  GEOMETRY 
FIGURE  4  a 


inches 


EDGE  INJECTION  GEOMETRY 
FIGURE  4b 
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DESCRIPTION  OF  THE  COMPUTER  CODE 


the  computer  code  used  to  perform  the  parametric  studies  is 


termed  NASOIN  which  is 
INjection.  The  equations 
dimensional  axisymmetric 
equations.  These  equations 


3U  .  3F  .  1  3 ( rG)  .nH 

TT  +  +  “n  1  _  . 


at 


ax 


3r 


an  acronym  for  Navier  Stokes  Jet 
employed  in  the  code  are  the  two 
unsteady  two  specie  Navier  Stokes 
are  given  below: 


n=0  Two  Dimensional  Flow 
n=l  Axi-symmetric  Flow 


pu  J 

pv 

pe 

pf 

G  - 

PU 

PU2 

+  °xx  ] 

puv 

+  Txr  . 

►  H  * 

( pe 

+  >*)U  +Txrv  +  qx  ' 

puf 

+  mx  1 

pv 

puv 

+ 

T 

pv2 

xr 

a 

( pe 

+ 

rr 

0 

pvf 

4- 

rr 

• 

m 

y 

)  v  +  T  u  +  q 
xr 


-o 


0 

0 

68 

0 

0 


where 


= 

den 

u 

» 

ax  i 

V 

= 

ve  r 

e 

=s 

in  t 

f 

e: 

spe 

m 

= 

mas 

r 

**> 

fraction 


stress  components 


This  code  has  evolved  over  a  number  of  years  and  is  based  on 
work  performed  in  references  4,5  and  6.  The  equations  are  solved 
using  the  explicit  technique  of  MacCormack  (7).  There  are  five 
independent  variables  which  are  solved  for,  the  density,  x 
velocity,  y  velocity,  total  energy  and  the  specie  mass  fraction. 
The  flow  is  not  considered  to  have  a  constant  total  temperature. 
Therefore  at  the  inflow  plane  or  base  injection  plane  the  static 
temperature  is  specified  as  a  function  of  y.  Inis  results  in 
a  total  temperature  and  therefore  a  total  energy  variation  in  the 
vertical  direction. 


CODE  MODIFICATIONS 


'  J  W  FV  U  ,*4  '<■«’»«  r  1 


m*n.  JC"'  ■*./•.  n^r 


A  significant  modi f icat ion  that  was  required  to  the  code  was 
implementing  a  subsonic  inflow  boundary  condition  for  the 
injected  mass.  This  was  done  using  an  additional  subroutine.  Two 
approaches  were  tried,  the  first  was  to  assume  a  total  injection 
pressure  that  was  only  clightly  larger  than  the  static  pressure 
in  the  base  region.  Using  these  two  pressures  the  injectant  Mach 
number  was  then  computed.  The  injectant  mass  flow  ratio  "I"  was 
specified  and  also  the  injectant  total  temperature.  Using  the 
total  temperature  and  the  Mach  number  the  injectant  static 
temperature  and  velocity  was  then  compjted.  Using  the  static 
pressure  and  the  static  temperature  the  injectant  density  was 
then  calculated.  An  estimate  for  the  mass  injectant  flow  rate 
ratio  "Ip"  was  then  computed  using  the  velocity,  density  and 
injectant  area.  This  value  was  then  compared  to  the  specfied 
injectant  mass  flow  rate  value.  The  first  value  was  always  less 
than  the  specfied  value  (the  initial  value  for  the  injectant 
total  pressure  was  chosen  so  that  this  was  always  the  case).  The 
injectant  pressure  was  then  increased  gradually  until  the 
calculated  value  for  the  mass  injection  was  greater  than  the 
specified  value.  The  correct  value  for  the  injectant  conditions 
was  then  interpolated  on.  These  values  were  then  used  to  compute 
a  new  estimate  for  the  mass  injectant  parameter.  If  the  two 
values  agreed  to  within  2%,  the  iteration  was  stopped,  if  not  the 
total  pressure  increment  for  the  injectant  was  decreased  and  the 
process  was  started  over  again.  This  approach  appeared  to  work 
until  at  about  9000  iterations  it  broke  down.  Large  values  of 
base  pressure  were  obtained.  The  reason  for  the  occurence  was  not 
resolved.  Because  of  this  a  second  approach  was  also  pursued. 

The  second  approach  was  to  express  the  mass  flow  as  a 
function  of  the  injectant  conditions  using  the  expression  for  the 
injection  mass  flow  ratio. 


I  =  lP.u.AJ; 

P-U-Ab 

Using  the  equation  for  a  perfect  gas,  the  definition  of  the  Mach 
number,  and  speed  of  sound  an  expression  for  the  injectant  Mach 
number  can  be  obtained. 


M  . 
J 


I 


A 


b 


Note:  The  pressure  used  in  this  equation  was  an  average  pressure 
at  the  first  flow  field  interior  grid  point  over  the  injectant 
exit  area. 
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This  expression  is  solved  for  the  Mach  Number  using  an  iteration 
process.  A  value  for  the  injectant  total  temperature  is  specified 
and  assumed  to  be  constant.  An  estimate  for  the  injectant  static 
temperature  (Tj  =  T0j )  is  assumed.  The  above  expression  is  used 
to  compute  the  Mach  number.  The  static  temperature  is  updated 
using  the  new  value  of  Mach  Number.  A  second  value  for  the  Mach 
Number  is  determined  and  compared  to  the  first.  This  process  is 
repeated  until  there  is  no  change  in  the  value  for  the  Mach 
Number.  Typically  this  required  8  iterations.  After  the  Mach 
Number  is  determined  the  static  temperature  is  computed,  the 
injectant  velocity  and  then  the  injectant  density.  This  second 
approach  was  successful. 

The  second  routine  that  was  added  to  the  code  was  one  that 
computes  the  force  on  the  projectile  base.  The  integration  of  the 
base  pressure  used  a  simple  trapezoidel  integration  technique. 
Separate  values  for  the  static  force  and  injectant  momentum  were 
computed.  The  injectant  momentum  force  is  not  include  in  the  base 
force  results  and  would  contribute  an  additional  10%  to  the  base 
force  for  injectant  values  greater  than  I=.04. 

In  addition  to  the  above  two  modifications  to  the  computer 
code  a  major  restructur ing  of  the  code  was  performed.  This  major 
restructur ing  was  to  replace  all  the  double  array  subscripted 
variables  with  a  single  array.  The  effect  of  this  is  to  increase 
the  vector  length.  During  the  process  of  performing  the  code 
restructur ing  any  coding  that  inhibited  do  loop  vectorization  was 
removed  from  the  loop. 

The  major  differences  between  the  two  versions  of  the  code 
are  that  the  single  array  version  has  no  'IF'  statements  inside 
the  inner  DO  loops  and  has  a  vector  length  of  NNX  (number  of 
points  in  the  x  direction)  times  NNY  (number  of  points  in  the  y 
direction).  For  example  for  a  100  (NNX)  by  50  (NNY)  grid  the  old 
version  had  a  maximum  vector  length  of  100  (the  x  dimension 
length).  The  new  version  now  has  a  vector  length  of  1500 
( NNX*NNY ) .  This  allows  use  of  the  maximum  length  vector  permitted 
by  the  CRAY  processors. 

The  changes  to  the  code  where  made  on  the  IBM  AT.  The  code 
was  then  was  compiled  on  the  IBM  PC  AT  to  locate  and  correct  any 
FORTRAN  syntax  errors .  Typically  it  took  about  20  minutes  to 
compile  the  4500  lines  of  code  on  the  IBM  AT.  The  code  was  then 
loaded  onto  the  CDC  CYBERNET  system  and  benchmarks  between  the 
old  version  and  the  new  single  array  version  were  made.  The 
results  from  the  single  array  version  were  compared  to  results 
from  a  previous  version  of  the  code.  Both  sets  of  results  were  in 
agreement  indicating  that  there  were  no  coding  errors  introduced 
as  a  result  of  the  code  changes. 


The  results  of  these  benchmarks  are  given  in  Table  2. 
Typically  the  new  single  array  version  runs  about  three  times 
faster  on  the  Cray  XM/P  and  three  times  faster  on  the  CYBER  205. 
This  new  version  of  the  code  does  not  have  the  subsonic  injection 
boundary  condition  routine  or  the  base  force  routine. 


Table  2 

Benchmark  Comparison 
(Performed  on  CDC  Cybernet  System) 


CPU  Execution  Time  in  Seconds 


Version 

Or  igina 1 

Single 

Optimized  Single 

Computer 

Code 

Array 

Array 

CYBER  205 

337.6 

110.1 

— 

CRAY  X-MP/24 

80-90* 

46.85 

28.6 

(These  times  do  not  include  compilation  time) 


CPU  Breakdown  on  Cray  X-MP/24 


Routine 

Single 

Optimized  Single 

Array 

Array 

NAS JIN 

1.246 

1.209 

TRAN 

.120 

.116 

THERMO 

4.667 

.964 

BC 

2.690 

2.634 

SIDE 

.873 

.861 

STP 

16.560 

3.450 

SOL  VR 

8.030 

8.116 

SMOOTH 

4.044 

3.239 

STRESS 

5.446 

5.166 

FLUXES 

2.180 

1 .916 

PRINT 

.905 

.668 

Other 

.089 

.061 

Total  CPU  time  (seconds)  46.85  28.6 

Benchmark  case  was  a  60x40  grid  run  for  1000  time  steps 

*  Estimated:  The  original  version  required  139  seconds  on  a 
CRAY-1  S/2000.  The  CRAY  X-MP/24  was  not  avaiable  for  this 
original  benchmark. 


COMBUSTION  MODEL 


The  Phase  I  effort  considered  the  effects  of  combustion  by 
employing  a  flame  sheet  hydrogen  air  combustion  model  (8).  In 
this  type  of  combustion  model  the  chemical  kinetics  employed  are 
those  of  chemical  equilibrium  or  an  infinitely  fast  chemical 
reaction.  This  is  a  useful  model  of  a  single  constituent  fuel 
combustion  process  that  is  essentially  a  simplified  treatment  of 
local  chemical  equilibrium.  This  identical  approach  has  been 
employed  in  references  5  and  6  for  analysis  of  a  hydrogen  fueled 
scramjet  combustor.  The  model  can  be  illustrated  by  considering  a 
hydrogen  air  system.  Restricting  the  treatment  to  problems  with  a 
maximum  temperature  of  less  than  2500°K,  an  examination  of  the 
pertinent  equilibrium  constants  reveals  that  it  is  reasonable  to 
neglect  all  reactions  involving  nitrogen  and  that  the  existence 
of  the  radicals  of  oxygen  and  hydrogen  (i.e.,  0,  H  and  OH)  can  be 
neglected.  Thus,  only  the  simple  overall  reaction 

H2  +  1/2  02  *  H20 

with  the  related  equilibrium  constant 


must  be  considered.  Furthermore,  for  T  <  2500°  K,  Kp,®  »  1 ,  so 
it  can  be  asserted  that  either  the  concentration  of  hydrogen  or 
oxygen  must  be  essentially  zero  in  certain  regions  of  the  flow. 
Thus,  we  come  to  the  flame  sheet  model  where  the  flow  is  divided 
into  two  regions:  one  where  there  is  no  fuel  and  one  where  there 
is  no  oxygen.  The  boundary  between  the  two  is  the  "flame  sheet" 
where  the  concentration  of  both  hydrogen  and  oxygen  is  zero.  This 
"flame  sheet"  occurs  at  the  locus  of  points  where  the  ratio  of 
oxygen  atoms  to  hydrogen  atoms  is  stoichiometric.  This  model  is 
illustrated  in  Figure  5. 

The  flame  sheet  model  is  an  equilibrium  chemistry  model. 
Therefore,  a  flame  sheet  dynamic  equilibrium  combustion  model  can 
be  considered  to  be  an  extreme  case  where  the  highest  degree  of 
heat  release  is  obtained  both  from  a  chemical  kinetics  and  fluid 
dynamics  model.  This  would  then  represent  one  extreme  for  the 
flow  phenomena  under  investigation. 


FLAME  SHEET  SCHEMATIC 


FIGURE  5 


The  method  of  solution  for  the  flame  sheet  model  model  is  to 
first  compute  the  flow  field  fluid  dynamics  for  a  single  time 
step.  This  determines  the  amount  of  fuel  that  is  diffused  and 
convected  by  the  flow  field  during  this  time  step.  Next  the 
amount  of  oxygen  and  nitrogen  present  at  each  grid  point  is 
determined.  This  is  done  using  the  relations: 


=  (1.0-  Vh2  )  *  .232  and 


=  1.0-YMa:-Yc 


The  .232  value  in  the  above  relation  is  the  gram  atom  weight  of 
oxygen  present  in  the  amount  of  air  and  is  determined  by  the 
following  calculation: 

Gram  atom  weight  of  air  =  W* 

Wa  =  Number  of  oxygen  atoms  x  oxygen  atomic  weight  + 

Number  of  nitrogen  atoms  x  nitrogen  atomic  weight  x 
number  of  nitrogen  atoms  for  each  atom  of  oxygen 

=  2  X  (16)  +  2  x  14.0  x  3.76  =  137.28 

Woa/Wo  =  32/137.20  =  .232 

The  local  stoichiometric  ratio  at  each  grid  point  is  determined 
based  on  the  local  hydrogen  and  oxygen  composition: 

Ys  =  (2.016/16.00)  Yoa 

A  comparison  is  made  to  determine  if  the  amount  of  hydrogen 
present  is  greater  or  less  than  Ys  .  If  the  amount  of  hydrogen 


present  is  greater  than  the  stoichiometric  value  the  amount  of 
hydrogen  consumed  in  the  reaction  is  set  equal  to  the 
stoich iometr ic  value.  The  residual  amount  of  hydrogen  is  then: 

(  Yh2  )■*  1  n»l  =  (  Vh2  )lnlfcl«l  ~  Yj 

If  the  amount  of  hydrogen  present  is  less  than  the  stoichiometric 
value  then  all  of  the  hydrogen  is  consumed  in  the  reaction  and 

(YHl)fln«l  =  0.0 

In  the  case  when  there  is  excess  hydrogen,  the  amount  of  oxygen 
and  water  are  computed  from  the  equation: 

( Yo2 )f  i n*l  =  (Yo  'linltlal  ~  1  6  .  0  * Y j  / 2 . 0 1 6 

(YH2o)fin.1  =  18.016  *  Ys/2.016 

or  for  the  case  when  there  is  not  enough  hydrogen  for  complete 
combust  ion 

(  Yo2  )  ■*■  i  n«l  =  (  Y02  ~  1  6 . 0  *  Yh2  /  2  -  0  1  6 

(  Yh2D  )  •»  1  n «  1  =  18.016  *  Ym2/2.016 

The  results  of  the  above  proceedure  allow  computing  the  gas 
composition  at  each  grid  point.  This  composition  is  then  used  to 
compute  the  static  enthalpy  of  the  mixture.  The  enthalpy  is 
then  used  in  the  equation  for  total  energy  to  compute  a  new 
value  of  temperature  after  the  fuel  and  air  has  burned. 


PARAMETRIC  CASES 

The  effects  that  were  investigated  consisted  of  determining 
the  impact  on  base  drag  of: 

1)  edge  injection  in  the  axial  direction 

2)  center  injection  in  the  axial  direction 

3)  edge  injection  at  an  angle  of  45° 

4)  effect  of  injectant  temperature 

5)  effect  of  spin 

6)  air  and  hydrogen  injection  with  no  reaction 

7)  air  and  hydrogen  injection  with  reaction 

The  effect  of  injectant  conditions  on  the  injection  mass 
flow  requirements  was  estimated  as  a  function  of  free  stream  Mach 
Number.  Results  of  these  calculations  are  shown  in  Figures  6,  7 
and  8.  These  curves  were  generated  only  to  obtain  a  relative 
comparison  of  different  conditions.  It  was  assumed  for  all  of 
these  cases  that  the  base  pressure  was  constant  and  equal  to  7.0 
psia  and  the  base  area  to  injectant  area  ratio  was  constant  at  a 
value  of  2.5.  The  figures  indicate  that  the  higher  the  injectant 
temperature  and  molecular  weight  the  lower  the  injectant  mass 
flow  requirements  to  obtain  the  same  value  of  base  pressure. 
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THEORETICAL  RESULTS 


EFFECT  OF  TEMPERATURE 
FIGURE  8 


The  majority  of  the  theoretical  results 
NAS0IN  code  were  performed  on  the  Ballistic 
CRAY  XMP/48.  The  computer  was  accessed  over  a 
line  at  1200  baud.  Typical  run  times  using 
of  the  code  required  30  minutes  to  90  minutes 
number  of  grid  points  and  the  number  of  time 
run  typically  was  for  20,000  time  steps  using 
points  in  the  x  direction  and  55  grid  points 
The  turn  around  time  was  either  overnight 
submitted  early  in  the  day  it  would  be 
afternoon . 


obtained  from  the 
Research  Laboratory 
toll  free  dial  in 
the  original  version 
depending  on  the 
steps.  The  90  minute 
a  grid  with  100 
in  the  y  direction, 
or  if  a  job  was 
finished  by  late 


The  theoretical  effort  performed  focused  on  determining  the 
effects  of  central  injection  and  edge  injection  on  projectile 
base  drag.  Prior  to  performing  the  injection  studies,  it  was 
desirable  to  correlate  the  results  of  the  code  with  experimental 
data  when  there  was  no  base  injection.  The  computed  base  pressure 
with  no  injection  are  plotted  in  Figure  9.  Also  shown  in  the 
figure  are  experimental  results  from  NACA  report  1051  "An 
Analysis  of  Base  Pressure  at  Supersonic  Velocities  and  Comparison 
with  Experiment"  by  Dean  R.  Chapman.  The  computed  values  of  base 
pressure  agree  to  within  10%  of  the  experimental  values.  The  base 
pressures  with  no  injection  from  AIAA  Paper  86-0487  "Supersonic 
Flow  over  Cylindrical  Afterbodies  with  Base  Bleed"  by  J.  Sahu  are 
also  shown  on  the  figure.  Again  these  values  agree  to  within  10% 
of  the  experimental  values.  The  original  data  from  NACA  1051  is 
shown  in  Figure  10  and  covers  a  wide  range  of  experimental  test 
conditions  and  body  configurations.  Therefore  agreement  to  within 
10%  is  not  unreasonable. 

Having  established  a  confidence  level  of  the  code  when  there 
is  no  injection,  the  objective  was  to  determine  the  confidence 
level  of  the  code  with  injection.  The  base  pressures  for  cases 
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with  injection  again  from  AIAA  Paper  86-0487  "Supersonic  Flow 
over  Cylindrical  Afterbodies  with  Base  Bleed"  are  shown  in  Figure 
11.  The  results  are  for  a  Mach  1.8  center  injection  case  whereas 
the  AIAA  results  are  for  a  Mach  nunber  of  1.7.  Although  the  body 
sizes  and  free  stream  conditions  are  not  the  same  for  these  two 
cases  there  is  general  agreement  between  the  two  sets  of  results. 
The  minor  difference  between  the  current  calculations  and  the 
previous  ones  is  that  there  is  no  minimum  value  for  base  drag  as 
the  injection  mass  ratio  "1“  increases.  The  base  drag  continues 
to  decrease  as  injection  increases,  indicating  the  more  mass 
injected  the  lower  the  base  drag. 
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Parametric  calculations  were  then  performed  for  center  and 
edge  injection  on  the  5.0  inch  diameter  body  of  Figure  4.  Results 
for  air  injection  are  shown  in  Figure  12  for  central  injection 
and  Figure  13  for  edge  injection.  Injectant  total  temperature  for 
these  cases  was  the  same  as  the  free  stream  total  temperature. 
Both  the  value  of  the  force  and  the  base  drag  coeffficinet  are 
plotted  in  the  figures.  In  general  for  center  injection,  as  the 
amount  of  mass  injection  increases  the  base  drag  decreases.  This 
is  true  for  both  central  and  edge  injection.  An  exception  to  this 
is  the  center  injection  Mach  1.4  case  which  exhibits  a  definite 
peak  at  I=.03.  The  edge  injection  results,  however,  indicate  that 
at  the  lower  Mach  number  of  1.4  a  substantial  drag  reduction 
occurs  with  only  a  small  amount  of  mass  injection  (I=.0l).  In 
fact  the  edge  injection  actually  results  in  forces  large  enough 
to  generate  some  degree  of  thrust.  As  the  free  stream  Mach  number 
increases  the  base  drag  coefficient  decreases  but  is  still  large 
enough  so  that  there  is  almost  no  base  drag  penalty  incurred.  In 
general  the  edge  injection  results  indicate  a  potential 
improvement  in  projectile  range  may  be  obtained  over  that  of 
central  injection.  To  determine  the  degree  of  improvement 
requires  performing  trajectory  simulations  for  the  different 
injection  techniques. 
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CENTER  INJECTION 
FIGURE  12 


A  possible  explanation  of  this  is  that  edge  injection  is 
more  effective  in  generating  an  equivalent  streamlined  body  then 
center  injection.  This  is  depicted  in  Figure  14a  and  14b.  For  the 
same  amount  of  mass  flow  the  center  jet  injectant  has  to  spread 
further  than  the  edge  jet  injectant  to  attenuate  the  expansion  of 
the  free  stream  gas  into  the  centerline  region  of  the  body. 
Therefore,  based  on  these  initial  results,  it  appears  that  edge 
injection  is  more  effective  in  the  reduction  of  base  drag  then 
central  injection.  These  results  however  should  be  verified  by 
additional  calculations  and  experimental  wind  tunnel  tests. 
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CENTER  INJECTION 
FIGURE  14a 
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EQUIVALENT  BODY  CONFIGURATIONS 
FIGURE  14 


The  effect  of  injection  angle  angle  for  the  edge  injection 
case  was  also  investigated.  Shown  in  Figure  15  are  the  results 
for  injecting  at  an  angle  of  +45°  and  -45°.  The  results  indicate 
that  angled  injection  in  general  detracts  from  the  reducing  the 
base  drag.  There  is  also  a  crossover  where  a  negative  injection 
angle  is  Initially  better  than  positive  injection  until  at  an 
Injection  ratio  of  I=.04  where  this  trend  is  reversed.  The 
reason  for  lookin  at  this  effect  was  based  on  the  fact  that  if 
the  injectant  was  directed  at  an  angle  of  45°  up  into  the 
oncoming  free  stream  flow,  it  would  appear  as  blockage  to  this 
incoming  stream  and  cause  an  increase  in  the  projectile  lip 
pressure  which  would  in  turn  cause  an  increase  in  the  base 
pressure.  To  fully  study  this  effect  a  three  dimensional  computer 
code  is  required  that  would  be  capable  of  analyzing  discreet 
injection.  The  area  between  the  discrete  injector  orifices  would 
allow  the  higher  projectile  lip  pressure  to  be  communicated  into 
the  base  region,  thereby  increasing  base  pressure  and  decreasing 
base  drag. 
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The  effect  of  injectant  temperature  was  then  investigated  at 
a  free  stream  Mach  number  of  1.8.  The  injectant  temperature  was 
increased  to  1100°  R  from  854°  R.  At  the  injectant  plane  (the 
base  region  station)  the  injectant  static  temperature  is 
specified  and  does  not  have  to  be  equal  to  the  free  stream  static 
temperature.  This  difference  in  static  temperature  at  the 
injectant  plane  also  gives  rise  to  a  variation  in  the  total 
energy  between  the  external  stream  and  the  injectant.  The  results 
of  these  calculations  are  shown  in  Figure  16.‘In  contrast  to  the 
previous  cases  where  the  injectant  total  temperature  was  equal  to 
the  free  stream  total  temperature,  there  is  now  a  more  distinct 
optimum  injectant  mass  flow  rate  that  occurs  at  I  =  .03.  The 
effect  of  increased  injectant  temperature  gives  a  further 
reduction  in  base  drag  than  cold  air  injection. 
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FIGURE  16 


The  effects  of  spin  on  central  injection  is  illustrated  in 
Figure  17  and  for  edge  injection  in  Figure  18.  The  spin  rates 
were  20,000  rpm  and  30,000  rpm  for  the  center  injection  case  and 
10,000  and  20,000  rpm  for  the  edge  injection  case.  Spin  causes  a 
lowering  of  the  base  pressure  as  is  shown  in  the  figures, 
apparently  the  spin  imparts  a  radial  outward  velocity  to  the  gas 
in  the  base  region.  This  decreases  the  density  of  the  base  region 
gas  and  therefore  also  the  base  pressure.  The  central  injection 
shows  a  monotonic  variation  of  the  base  force  with  spin.  The 
effect  of  spin  on  edge  injection  is  to  decrease  the  erratic 
behavior  that  occurs  with  no  spin. 
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EFFECT  OF  SPIN  CENTER  INJECTION 
FIGURE  17 


The  effects  of  injecting  a  mixture  of  hydrogen  and  air  for 
the  center  injection  case  is  shown  in  Figure  19  and  for  the  edge 
case  in  Figure  20.  The  amount  of  hydrogen  that  was  injected  is  a 
mass  fraction  of  .007.  Although  this  may  seem  small  the  hydrogen 
mass  fraction  for  a  stoichiometric  air  hydrogen  mixture  is  .027. 
Thus  the  amount  of  hydrogen  injected  is  about  one  fourth  the 
stiochiometr ic  value.  This  value  was  choosen  to  avoid  an 
excessive  amount  of  heat  release  that  may  cause  the  code  to  go 
unstable.  The  results  in  the  figures  indicate  that  the  hydrogen 
causes  the  base  force  to  increase  and  base  drag  to  decrease. 
This  agrees  with  other  published  work  which  indicates  that  the 
injection  of  light  molecular  weight  gases  in  the  base  region  are 
more  effective  in  decreasing  base  drag. 


INJECTION  MASS  RATIO  -  I 


CENTER  INJECTION  WITH  HYDROGEN 


FIGURE  19 


The  effects  of  combustion  in  the  base  region  is  shown  in 
Figure  21  for  the  edge  injection  case.  The  amount  of  hydrogen 
that  was  injected  is  the  same  as  in  the  previous  case  (hydrogen 
mass  fraction  =  .007).  The  results  in  the  figure  indicate  that 
the  combustion  of  hydrogen  in  the  base  region  caused  the  base 
force  to  decrease  and  therefore  the  base  drag  to  increase.  This 
is  unexpected  since  combustion  should  cause  an  increase  in 
temperature  and  any  increase  in  temperature  would  give  rise  to  a 
decrease  in  base  drag.  Possibly  the  amount  of  hydrogen  injected 
is  too  small  and  since  it  is  located  near  the  edge  of  the 
projectile  the  combusted  gas  mixes  rapidly  and  is  in  effect 
quenched  by  the  outer  cooler  stream. 


The  effects  of  combustion 
Figure  22  for  the  center 
hydrogen  that  was  injected  is 
The  results  in  the  figure 
injection,  the  combustion  of  a 
base  region  initially  causes 
and  therefore  the  base  drag  to 
flow  increases,  however,  thi 
injection  ratios  greater  than 
for  the  no  injection  case. 


in  the  base  region  is  shown  in 
injection  case.  Again  the  amount  of 
the  same  as  in  the  previous  cases 
indicate  that  for  the  central 
small  amount  of  hydrogen  in  the 
the  base  force  to  increase  (I=.01) 
decrease.  As  the  amount  of  mass 
s  trend  reverses  itself  and  at  mass 
.025  the  base  force  is  smaller  than 


In  obtaining  these  results  the  base  force  did  not  attain  a 
steady  state  value.  The  base  force  time  history  for  these  four 
injection  cases  are  shown  in  Figure  23.  All  of  the  cases 
demonstrated  an  almost  constant  value  of  base  force  at  early 
values  of  time.  At  later  values  of  time  this  force  increased.  It 
is  not  known  if  this  behaviour  is  due  to  improper  imp lemenat ion 
of  the  injectant  boundary  conditions,  downstream  boundary 
conditions  or  reflections  from  the  downstream  boundaries  that 
traveled  back  to  the  base  region.  Each  of  these  effects  can  be 
investigated  by  moving  the  downstream  boundary  or  using  different 
boundary  conditions  such  as  characteristics  rather  than 
extrapolated  boundary  conditions. 
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A  problem  area  that  surfaced  during  the  performance  of  the 
theoretical  calculations  was  the  base  force  in  some  cases  did  not 
readily  converge  to  a  constant  value.  The  degree  of  oscillation 
is  depicted  in  Figure  24  for  a  case  with  no  injection  and  in 
Figure  25  for  a  case  where  the  injectant  mass  flow  ratio  was  .01 
and  the  free  stream  Mach  number  was  2.2.  This  oscillation  is 
believed  to  be  due  to  applying  a  constant  value  of  pressure  over 
the  injection  plane  where  in  reality  the  pressure  and  velocity 
vary  over  the  injection  plane.  This  effect  or  variation  would  be 
more  dominant  for  large  bodies  with  sparser  grids  than  for 
smaller  diameter  bodies  with  more  tightly  clustered  grids. 
Comparison  to  the  no  injection  case  shows  that  the  base  force 
converges  very  rapidly  to  a  constant  value  after  only  about  4000 
iterations.  Except  for  the  peak  that  occurs  at  about  5000 
iterations  the  value  is  extremely  steady.  This  indicates  that 
further  work  must  be  done  to  either  smooth  the  flow  in  the  base 
injection  region  or  to  allow  for  a  variable  injectant  conditions 
at  the  injectant  plane. 


BASE  FORCE  TIME  HISTORY 
FIGURE  24 


Conclusions 


The  results  obtained  indicate  that  subsonic  base  injection 
can  be  beneficial  in  reducing  base  drag.  The  use  of  edge 
injection  gives  higher  values  of  base  drag  reduction  (by 
approximately  20%)  than  center  injection.  Use  of  higher 
temperature  injectant  gas  also  gives  larger  values  of  base  drag 
reduction.  It  appears  that  small  injection  amounts  of  burning  gas 
in  the  base  region  is  more  beneficial  than  large  amounts,  which 
might  blow  off  the  base  bubble. 

Details  of  the  optimum  injection  scheme,  however,  still 
remain  to  be  determined.  In  particular,  effects  of  flow  channels 
located  in  the  base  region  that  direct  the  flow  either  radially 
inward  or  outward  may  increase  or  possibly  decrease  the  base 
pressure  and  therefore  alter  the  range  of  the  artillery  shell. 
Potential  flow  channel  concept  are  illustrated  in  Figure  26.  The 
objective  of  these  concepts  are  to  direct  injectant  gas  so  that 
it  interacts  with  ambient  air  that  is  trying  to  turn  the  corner 
or  to  use  bleed  paths  to  generate  small  jets  which  interact  with 
the  gaseous  injectant. 
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mass  injection  ratio 
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viscos i ty 
dens i ty 

angular  velocity 


SUBSCRIPTS 


base  value 
injectant  value 
injectant  supply  value 
free  stream  condition 


43 


